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Exact asymptotics of the current in boundary dissipated quantum chains in large external fields 
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Boundary driven quantum master equation for a general inhomogeneous (non-integrable) anisotropic Heisen¬ 
berg spin 1/2 chain, or an equivalent nearest neighbor interacting spinless fermion chain, is considered in the 
presence of a strong external field /. We present an exact closed form expression for large / asymptotics of 
the current in the presence of pure incoherent source and sink dissipation at the boundaries. In application we 
demonstrate arbitrary large current rectification in the presence of interaction. 
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Introduction - The study of quantum dynamics far from 
equilibrium is currently at the forefront of theoretical m 
(but also experimental EHU) condensed matter physics re¬ 
search. Although most of non-equilibrium phenomena stud¬ 
ied are directly time dependent, such as e.g. in the so-called 
quantum quenches 0), there is a growing interest in study¬ 
ing steady-state phenomena eoed especially under large ther¬ 
mal/chemical bias and/or external fields BED- A prominent 
examples of this sort are dielectric breakdown in the Mott 
insulator HMD or rectification effects in nanoscale heat 
or particle transport fT9H2TI . Particularly fruitful direction 
which has been developed recently ll22ll is a hybrid approach 
of open quantum systems formalism l23l . where the bulk of 
a finite (possibly large) interacting quantum system is treated 
fully coherently with a unitary generator, whereas the environ¬ 
ment is considered to act locally on the system’s boundaries in 
terms of incoherent markovian quantum noise. Remarkably, it 
has been demonstrated that such an approach allows for inter¬ 
esting exactly solvable, integrable instances and can 

be phenomenologically justified even for non-small system- 
environment couplings in terms of the so-called repeated in¬ 
teractions protocol l28l . 

An exact treatment of nonequilibrium transport has been 
so far possible only for pseudo-force driving originated from 
the bias in markovian processes at the system’s boundaries. 
In this Letter we combine this approach with a real force f, 
resulting from an external field gradient, and formulate a sys¬ 
tematic perturbation theory of the steady state current in 1//. 
We show that under a plausible assumption about the zeroth 
order of the steady state density matrix an explicit asymp¬ 
totic f —y oo expression for the current can be derived for 
a general inhomogeneous (and nonintegrable) XXZ spin 1/2 
chain. For long chains one can demonstrate arbitrarily large 
relative current rectification under reversal / —» —/. Our 
formula is clearly corroborated by numerical computations 
on small chains, where further resonance effects of unusually 
large currents at particular values of / are discussed and partly 
explained. The result applies as well for large anisotropy (in¬ 
teraction) asymptotics. 

Quantum transport in boundary & field driven spin chain.— 
We consider an open inhomogeneous XXZ spin 1/2 chain 
exposed to an additional magnetic field, so that bulk of the 


system is governed by the Hamiltonian H = H xxz + V, 

n—1 

H X xz = H xx -\- H Z) H xx = 'y ^ e loLj{crj (jj _)_]_ (Jj 

3 =1 

n—1 n 

^ = Z A i*M + 1, v = Y,frt- (1) 

7=1 7=1 

where parameters oij , A j may be site dependent. For the ho¬ 
mogeneous case ctj = 1, A ? = A will be used, fj determine 
the profile of magnetic field. Using Wigner-Jordan transfor¬ 
mation the problem can be mapped to spinless fermion chain 
where otj map to free hopping amplitudes, A ; to interactions 
and fj to electric potential. Density matrix p of the system 
may be described using the Liouville master equation in the 
Lindblad form |22] [23j 

= Cp := —i[H, p{t)} + t>p, (2) 

where Vp = 2 L^p(t)L\ — {L^L^, p(t)} describes cou¬ 

pling to the bath at the edges. Here, only the extreme case of 
pure source/sink on the left/right end, L\ = \fza'\, 1 .2 = 
y/ea~, will be considered. Nonequilibrium steady state 
(NESS) poo = Hindoo p(t) is found as the fixed point 
Cpoo = 0. The effect of driving will be characterized by the 
spin current from site j to j + 1 

Jj = iaj ( 0 + a~ +1 - erf a + +1 ), (3) 

which is in the limit t —► 00 , i.e. for NESS, site independent, 

J = tr{JjPoo)/ trpoo- 

Large field, large interaction asymptotics.- We split the Li- 
ouvillian into diagonal and off-diagonal parts C = fCo + £ 1 , 

Cop = -i[V + H z ,p\, £\p = -i[H xx ,p\ +Vp, ( 4 ) 

where H z = E"=i' A ja]a] +1 , V = £"=1 j/crj with A, = 
A j If, fj = fj //. We consider 1// as a formal small param¬ 
eter, assuming / aj, e. Note that both diagonal generators, 
V or //-, need not be large. We may either fix f n — fi = 1, 
in which case / represents the total field/potential drop across 
the system while A j = 0(1), or consider the large anisotropy 
(interaction) asymptotics A ? = 0(1) while fj = 0(f~ 1 ). 
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Expanding NESS in powers of / \ p^ = Y)^=o f p P < ' p \ 
different orders are connected with the recurrence relation: 

£oP (p+1) + Ap (p) = o, £ 0 p (0) = o. (5) 

Each order of p^ can be expanded in terms of a complete 
Pauli basis p (p) = E{ Sl ,..., Sn } II j where Sj e 

{0, z, +, —}, and a 0 = 1. Condition for p (°) in Eq. (J 5 J is sat¬ 
isfied for exponentially many products . a^ for which the 

sum E"=i fjtsj + E”=i A j ((&> 3 - + &,-+ 1 ) ( mod 2 )) van¬ 
ishes, where = ±1 for Sj = ± and £ S;j = 0 other¬ 
wise. However, exact numerical diagonalization of £ on small 
systems of even size n < 8 shows that for strong field or 
anisotropy p ! ' {y> is a diagonal projector 


O (o) 


j = 1 i=f+t 


with step-like magnetization profile. In case of large 
anisotropy this has already been observed |29| and given 
physical motivation suitable for our limit as well. Thus, we 
conjecture that p y}) has the form of Eq. (JfiJ for a system of 
arbitrary even length n. If n is odd p (°) contains two terms 
of the form 0 with a jump in magnetization between sites 
(" 2 ‘ 1 , ^UT) or ("+ 1 . "+ 3 ). Further on we assume even n, 
generalizing the solution to arbitrary n at the end of deriva¬ 
tion. 

We introduce the following notation for local operators 


i(i+»■) = (!). «- = ({)• 

5 ( 1 -o Mi), » + = (!). (7) 

i.e., (/) = \v)(p\, v,p£ {t4}. E.g., p is represented as 
a single component 


P {0) = ao 


t t 
t - f 


4 . ... 1 \ 
4 . ... 4 .)■ 


( 8 ) 


Vertical bar marks the center of the system. It is instructive 
to think of the effect that different parts of £ have on a cho¬ 
sen term MjIX; a j J ' n suc h notation. The action of £0 
is diagonal, i.e. it does not change the operator product to 
which it is applied, while unitary part of C \ causes exchange 
of neighboring pairs if different, (tjI. 7 + 1 ) °r )■ Dissi- 

pator T> is the only process to relate imaginary and real part of 
coefficients involved in the recurrence relation and is trivial if 
both p, v are aligned with the boundary driving. 

Using the recurrence relation { 5 } higher order p (p> can be 
generated from p(°\ E.g. the first order /Z 1 * then consists of 
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(9) 

However, we will not strive to built the whole NESS density 
matrix, instead we will relate perturbatively only its compo¬ 
nent relevant for the NESS spin current. 


Steady state spin current.- We would like to identify the 
leading order terms in the current J = tr (JjPoo)/ trpoo. In 
our diagrammatic notation tr is a sum of coefficients at con¬ 
figurations with the same upper and lower index at each site. 
Therefore trpoo = ao + 0(f~ 2 ). Similarly, trJjPaa is a 
sum over configurations in p^ with (t±<jJ +1 and 1; ±of else¬ 
where. Suppose we calculate current from the operator J n / 2 . 
We will show that then the leading contribution comes from 
the components with 1; -f <rf for l < n/2 and 1; — erf for 
l nf 2 ~b 1. 

We claim that in order to obtain the current from tr J ra / 2 Poo 
up to £>(/ -(n+2) ) only the terms p = p(°) + pi.R + pi^ + 
P 2 ,r + P 2 ,l in the density matrix are relevant. Expression 


Pi,R = 
a 1 

2/ f 



t 
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t 4- 11 4- 4- 4- j 
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with a p G 1R introduces a sequence of terms connected by the 
recurrence relation ([5j>, which takes upper tn/2+1 in P^ (in 
red color) to the right boundary and back to its initial position. 
After each step coefficient involved is for a / _1 higher in ex¬ 
pansion, being real until the action of V and imaginary after. 
pi t L is constructed symmetrically, taking upper i n i 2 to the 
left boundary and will be expanded in terms of coefficients b p 
with 61 = a\. Other p 2 .R(L) originate from the second term of 
/c/T in Eq. (pk and correspond to shifting lower tn/2+1 (in/ 2) 
to the right (left) boundary. 

From this construction it is clear that the leading contribu¬ 
tion to the current is of order f~ n , namely 


J= — 


v 77 2 [a, n + bn) 
f n a 0 


o(f 


~{n+ 2) 


( 10 ) 


since real coefficients a .\. h \ get subtracted, while nth order 
coefficient in pin and p 2 ,R (equivalently for p\i, and P 2 ,l) 
are complex conjugated and sum up yielding a factor 2 in 


Eq. (10 1 . Any other operator product term having o’„/ 2 cr n/ 2 +i 
and projectors 1* + < 7 ? for some i > n/2 +1 is of higher order 
in / -1 due to obvious additional steps needed to create it from 
p(°\ To calculate the spin current we have to express a n , b n 
with a ( |. In the procedure each step yields an equation that 
relates coefficients from the consecutive orders via recurrence 
0- Explicit relations obtained from pi R are: 
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Step 1: 


— 2ia:±a 0 + 2i + A~ +1 + /« — ft±+i) «i — 0 


a i = 


a 0 a; 


Af-i + A f+i ^ /f+i + /f 

Step A:, 1 < k < n/2: 

«fe -1 ctf+fc-t 


( 11 ) 


ak = 


An+fe-l + A™_|_fc + A»_l — A" — /fi>+fc + / 


Step n/2: 


Q , n_l Q? n _i 


A n _l + A»_l — A« — /n + 


( 12 ) 


(13) 


Step n/2 + 1 (action of dissipator V): 


am e 


a^+i = 


2 [A n _i + A»_i -A» -/n + /f] 


(14) 


Step n/2 + 1 + fc, l</c< n/2: 


«=+i+fc — 


^£+/c &n—k 


^n—k—1 “I - A n—k “t“ ^^ — 1 ^^ fn—k H” f V/ 

(15) 


Combining Eqs. EH T5j ) we finally obtain a n as a function 
of ao and parameters of Hamiltonian. b n is expressed analo¬ 
gously. 


— 

£ot n — 1 

Zk — l 

Vr Q f+fc-i “f+fc 

(16) 

a 0 

2(^A n -l+An _ ± — An —fn+f ™ 

\ 2 J. I / ~ ~ ~ ~ ~ ~ \2i 

) k = 1 { A % + k+ A % + k~l+An- 1 -An-fn +k + fn J 

— 

£OL\ 

-—1 

2 __ Ctrl _u _11 an _ k 

(17) 

do 

2(Ai+An +1 -Aa -/A+i-b/i) 

11 / ~ ~ ~ ~ ~ ~ \2 • 

k=1 (^Af-fc+Af-fc+l+A« + 1 -An -fn +1 +fn_ k+1 ) 


For a homogeneous chain and field which is odd with respect 
to the center /, = —f n +i~j expression for the current simpli¬ 
fies 

o f _1 

J = fn ( A - /» + fi) II ( 2A ~ 4+* + 4) ■ 

^ fc=l 

(18) 


One should note that a shift of fj for a constant does not 
change the result. Obtained expression for the current is ap¬ 
propriate for a general monotonic magnetic field, however, in 
the case of linear field profile fj — fj -1 = g it has especially 
simple form in terms of a gradient g 


T 8s ( 2A 
J = —- n 


9" 


9 


-2 7 


¥-1 


n 

k=1 


9 


-2 


(19) 


For systems of odd size n one should consider J( ra _i )/2 (or 
J( n+ i)/ 2 ). The leading order contributions originate from 
the part of p(°) with jump in magnetization between sites 
(t/A, ( or (tt±l. 2±3))_ The imaginary coefficient rel¬ 
evant for the current is then obtained as in system of even 
size 7 i—l, having contribution only from the boundary that is 
closer to the magnetization jump. The current for odd n, ho¬ 
mogeneous XXZ chain and field of odd symmetry w.r.t. the 
center of the system is then simply J(n) = J{n — l)/4. For 
an inhomogeneous chain or non-symmetric monotonic field 


the two components in p(°) have different weights, causing an 
additional prefactor in the latter equation. 

In Fig. [] we compare numerical and analytical results 
for the spin current J at linear field profile and homoge¬ 
neous Hamiltonian with A = 5 in systems of size n = 
4, 6 , 8 . At n = 4 we show also an example with in¬ 
homogeneous Hamiltonian with parameters {ai, 0 : 2 , 0 : 3 } = 
{0.1,0.3,0.5}, {Ai, A 2 , A 3 } = {-2.5,-5,-1}. Negative 
Aj are chosen for clarity since J{—Aj,—fj) = J(Aj.fj). 
Evidently, expression ( p~9| ) perfectly captures the dependence 
J (g) at large enough gradients g and indicates the position of 
some of the resonances in the numerical result. From the con¬ 
struction of the perturbative solution it is clear that poles in 
Eqs. ( | 1 6| 17 > appear when at some level of perturbative hier¬ 
archy the action of diagonal £ 0 °n considered operator prod¬ 
uct is trivial. Consequently, the hierarchy is disconnected and 
our expression not valid any more. Additional resonances ob¬ 
served in numerical results originate from other operator prod¬ 
ucts which are annulled by £ ( , that are not contained in p. 
In the resonances the hierarchy can be altered in a way that 
yields a contribution to the current in lower order of / -1 ||30l . 
The range of validity of asymptotic expression ( p~9| ), which 
is g £ (—00, 0) U (< 7 C , 00) for A > 0, is set by the furthest 
resonance, yielding a critical gradient g c = 2A. 

Flux rectification.- Asymmetry of the driving is revealed by 
changing the direction of the driving force, /—>■—/, result¬ 
ing in J{f) ^ J(—/). Rectifying behavior is conveniently 
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Figure 1. (Color online) Spin current J as a function of gradient g 
for linear field profile and anisotropy A = 5 and e = 1 in systems 
of size n = 4, 6, 8 obtained numerically (num) and from Eq. 03 
(ana). Curve at n = 4 labeled by A i„h shows same comparison for 
an inhomogeneous chain with parameters given in the text. 



Figure 2. (Color online) Dependence of flux rectification R on sys¬ 
tem size n for A > 0,e = 1 and gradients g = 3A,5A, 10A 
obtained from Eqs. ( |20|19[ >. Circles correspond to exact numerics 
for small systems. 


interaction prefers aligned spins, magnetic field with g > 0 
dictates the same orientation of spins as the Lindblad driv¬ 
ing, resulting in an enhanced insulating behavior of NESS. If 
g < 0 their effects are contrary, manifested in the enlarged 
spin current. Relation J(—Ay, —fj) = J(Aj, fj) generalizes 
the argument to A y > 0. 

For gradients above the critical value the rectification co¬ 
efficient R can be calculated for any system size n using the 
current asymptotics ( [T6p9| - In Fig.[2]we plot the dependence 
of R on the system size n (even only) for a homogeneous 
chain and linear field profile at three different gradients, g = 
3A,5A, 10A. Circles correspond to exact results for small 
systems. Our results clearly show that in the thermodynamic 
limit R reaches the extreme value R = sgn(A), cor¬ 

responding to insulating behavior for gradients of opposite 
sign than anisotropy A. 


Discussion - We have developed a systematic asymptotic 
expansion of the spin/particle current in quantum spin/fermion 
chains with nearest neighbor interaction driven with inho¬ 
mogeneous external fields and dissipated with incoherent 
source/sink baths at the chain ends. Our results have been 
worked out in detail on an example of an inhomogeneous 
XXZ spin 1/2 chain, driven by a transverse magnetic field 
gradient. We can point out two mechanisms for enhanc¬ 
ing current rectification: (i) rectification coefficient R —> 1 
in thermodynamic limit n —» oo, at fixed large gradient 
f/(n— 1), for any nonzero anisotropy/interaction A / 0, 
as a result of our asymptotic current formula; (ii) for reso¬ 
nant — field (gradient) values where our asymptotic analysis 
fails, the current is of higher order in 1 // for one direction of 
the field / than in the other. This clearly yields rectification 
coefficient R — > 1 even at fixed finite size n. Note that there 
is no rectification effect in noninteracting chains A = 0, in 
agreement with BTI . Our analysis also applies to reveal large 
anisotropy asymptotics at v anishin g external fields by setting 
Ay = /Ay, fj = 0 in Eqs. (16 17 1 , yielding for the homoge¬ 
neous XXZ chain J = e 2 _ 7 A _n , 7 = n 2 /4 — n /2 — 1 for 
even n and J = e 2 _< 5 A~( rl_1 ), 5= [n— l) 2 /4— (n— l)/2+l 
for odd n, as noted numerically in 129) . 

Discussions with M. Znidaric and E. Ilievski, as well as 
support by the grants PI-0044, J1-5349 of the Slovenian Re¬ 
search Agency are acknowledged. 


quantified by rectification coefficient lUTI 


( 20 ) 


which is zero in case of no rectification and reaches extreme 
values ±1 when the system for fields in one direction behaves 
as a perfect insulator. 

For one direction of the dissipative Lindblad driving the 
anisotropy Ay and the field /y have a cooperative effect and 
competing one for the other direction l32l . Fixing dissipative 
driving to our initial choice and assuming Ay < 0, when the 
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